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ABSTRACT 

Cosmological models that include suppression of the power spectrum of density fluctuations on small scales 
exhibit an exponential reduction of high-redshift, non-linear structures, including a reduction in the rate of gamma 
ray bursts (GRBs). Here we quantify the constraints that the detection of distant GRBs would place on structure 
formation models with reduced small-scale power. We compute the number of GRBs that could be detectable 
by the Swift satellite at high redshifts (z ^ 6), assuming that the GRBs trace the cosmic star formation history, 
which itself traces the formation of non-linear structures. We calibrate simple models of the intrinsic luminosity 
function of the bursts to the number and flux distribution of GRBs observed by the Burst And Transient Source 
Experiment (BATSE). We find that a discovery of high-z GRBs would imply strong constraints on models with 
reduced small-scale power. For example, a single GRB at z ^ 10, or 10 GRBs at z <; 5, discovered by Swift during 
its scheduled two-year mission, would rule out an exponential suppression of the power spectrum on scales below 
R c = 0.09 Mpc (exemplified by warm dark matter models with a particle mass of m x = 2 keV). Models with a 
less sharp suppression of small-scale power, such as those with a red tilt or a running scalar index, n s , are more 
difficult to constrain, because they are more degenerate with an increase in the power spectrum normalization, erg, 
and with models in which star-formation is allowed in low-mass minihalos. We find that a tilt of 5n s as 0.1 is 
difficult to detect; however, an observed rate of 1 GRB/yr at z <; 12 would yield an upper limit on the running of 
the spectral index, a = dnjdlnk > -0.05. 

Subject headings: cosmology: theory - dark matter - early Universe - galaxies - gamma rays: bursts: formation 
- galaxies: high-redshift - large-scale structure of universe 



1. INTRODUCTION 

In the years leading up to the recent launch of the Swift satel- 
lite, 1 it has been increasingly recognized that distant gamma 
ray bursts (GRBs) offer a unique probe of the high redshift 
universe. In particular, GRBs are the brightest known electro- 
magnetic phenomena in the universe, and can be detected up to 
very high redshifts (e.g. Wijers et al. 1998; Lamb & Reichart 
2000; Ciardi & Loeb 2000), well beyond the redshift z » 6.5 of 
the most distant currently known quasars (Fan et al. 2003) and 
galaxies (Kodaira et al. 2003). 

There is increasing evidence that GRBs are associated with 
the collapse of short-lived, massive stars, including the associ- 
ation of bursts with star-forming regions (e.g. Fruchter et al. 
1999), a contribution of supernova light to the optical afterglow 
(e.g. Bloom et al. 1999; Garnavich et al. 2003), and most di- 
rectly, association with a supernova (Stanek et al. 2003; Hjorth 
et al. 2003). 

As a result, the redshift distribution of bursts should fol- 
low the mean cosmic star-formation rate (SFR). Several studies 
have computed the evolution of the expected GRB rate under 
this assumption, based on empirical models of the global SFR 
(Totani 1997, 1999; 1999; Wijers et al. 1998; Lamb & Re- 
ichart 2000; Ciardi & Loeb 2000). Recent determinations of 
the cosmic SFR out to redshift z ~ 5 (e.g. Bunker et al. 2004; 

'See http://swift.gsfc.nasa.gov 



Gabasch et al. 2004; Giavalisco et al. 2004) show that star- 
formation is already significant at the upper end of the mea- 
sured redshift range, with J> 10% of all stars forming prior to 
z = 5, which would result in a significant population of GRBs 
at these redshifts. Further associating star-formation with the 
formation rate of non-linear dark matter halos, and using the- 
oretical models based on the Press & Schechter (1974) formal- 
ism, Bromm & Loeb (2002) and Choudhury & Srianand (2002) 
have extrapolated the SFR and obtained the GRB rates expected 
at still higher redshifts. These studies have concluded that a 
significant fraction (exceeding several percent) of GRBs de- 
tected at Swift s sensitivity should originate at redshifts as high 
as z > 10. The spectra of bright optical/IR afterglows of such 
distant GRBs can then reveal absorption features by neutral hy- 
drogen in the intergalactic medium (IGM), and can serve as an 
especially clean probe of the reionization history of the universe 
(Miralda-Escude 1998; Lamb & Reichart 2000; Choudhury & 
Srianand 2002; Lamb & Haiman 2003; Barkana & Loeb 2004). 

In this paper, we investigate a different method to utilize dis- 
tant GRBs, and to glean information about early structure for- 
mation. The mere presence of a GRB at, say, z > 10 will indi- 
cate that non-linear dark matter (DM) structures already exist at 
this redshift: the stars that give birth to the GRBs must form out 
of gas that collected inside dense DM potential wells. Structure 
formation in a cold dark matter (CDM) dominated universe is 
"bottom-up", with low-mass halos condensing first. In the cur- 
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rent concordance cosmology, with densities in cold dark mat- 
ter (CDM) and dark energy of (0 M ,^ a) ~ (0.3,0.7) that have 
emerged from WMAP and other recent experiments (Spergel 
et al. 2003), DM halos with the masses of globular clusters, 
10 5 ~ 6 M Q condense from ~ 3cr peaks of the initial primordial 
density field as early as z ~ 25. It is natural to identify these 
condensations as the sites where the first astrophysical objects, 
including the first massive stars, were born. As a result, one 
expects to possibly find GRBs out to this limiting redshift, but 
not beyond. 

With a scale-invariant initial fluctuation power spectrum, the 
CDM model has been remarkably successful, and has matched 
many observed properties of large-scale structures in the uni- 
verse, and of the cosmic microwave background (CMB) radi- 
ation. However, the power spectrum on smaller scales, corre- 
sponding to DM halo masses of M <, 10 9 M Q , remains poorly 
tested. Recent observations suggest that the standard model 
predicts too much power on small scales: it predicts steep cusps 
at the centers of dark matter halos, whereas the rotation curves 
of dwarf galaxies suggest a flat core; it also predicts more small 
satellites than appear to be present in the Local Group (these 
and other problems with CDM on small scales are reviewed by, 
e.g. Sellwood & Kosowsky 2001 and Haiman, Barkana & Os- 
triker 2001). Although astrophysical explanations of these ob- 
servations are possible, much recent attention has been focused 
on solutions involving the properties of dark matter. Propos- 
als include self-interacting dark matter (Spergel & Steinhardt 
2000), a repulsive interaction to gravity (Goodman 2000; Pee- 
bles 2000), the quantum-mechanical wave properties of ultra- 
light dark matter particles (Hu et al. 2000), and a resurrection 
of warm dark matter (WDM) models (Bode et al. 2001). 

By design, a common feature of models that attempt to solve 
the apparent small-scale problems of CDM is the reduction of 
fluctuation power on small scales. In addition, we note that 
reduced small scale power is a direct consequence of a range 
of slow-roll inflationary models, which predict a red tilt of the 
power spectrum, n s = d In P(k) / dink < 1, and a running of the 
spectral index a = dn s /d\nk < (see, e.g., Kinney 2003 for a 
general discussion of power spectra predicted in different infla- 
tion models, and Kosowsky & Turner 1995 for a discussion of 
models with a running index). Interest in such models was re- 
cently re-kindled, as they appeared preferred by a combination 
of CMB anisotropy data from WMAP with fluctuation statistics 
in the Lyman a forest (Spergel et al. 2003; Peiris et al. 2003). 

The loss of small-scale power generically suppresses struc- 
ture formation most severely at the highest redshifts, where the 
number of self-gravitating objects is drastically reduced. In 
each model, there exists a redshift beyond which the number 
of GRBs (or any other object) is exponentially suppressed, and 
a detection of a GRB beyond this redshift can be used to con- 
strain such models. A similar constraint can be obtained from 
the observed reionization of the universe at high-redshifts. For 
example, Barkana, Haiman & Ostriker (2001; hereafter BHO) 
showed that in the case of WDM models, invoking a WDM par- 
ticle mass of m x ~ 1 keV (approximately the mass required to 
solve the problems listed above), the paucity of ionizing sources 
makes it difficult to account for the reionization of the uni- 
verse by redshift z ~ 6. Reionization as early as z ~ 17, as re- 
cently suggested by WMAP observations of CMB anisotropics 
(Spergel et al. 2003) would require extreme efficiencies for star 
formation and ionizing photon production for masses ofm x < > 2 
keV (see also Somerville, Bullock & Livio 2003, who reach 



similar conclusions). 

GRBs, if discovered at z > 6, have the potential to provide 
independent and stronger constraints. The purpose of this paper 
is to quantify the constraints that the detection of distant GRBs 
would place on structure formation models with reduced small- 
scale power. Throughout most of our calculations, we focus on 
WDM model as a proxy, but our results are valid for any theory 
which imposes a small-scale cut-off in the primordial power 
spectrum. 

The rest of the paper is organized as follows. In § 2, we 
briefly describe our Monte Carlo approach to modify the stan- 
dard Press-Schechter formalism, allowing us to compute halo 
mass functions in WDM models. In § 3, we describe our method 
to compute the GRB rates using the halo mass functions, and 
simple models for the intrinsic GRB luminosity function. In 
§ 4, we present the constraints that high redshift GRB detec- 
tions would place on the small-scale power spectrum. Finally, 
in § 5, we discuss the implications of this work and offer our 
conclusions. 

Unless stated otherwise, throughout this paper we assume 
standard cosmological parameters, (Oa, fi M , n, erg, Ho) = 
(0.73, 0.27, 0.044, 1, 0.85, 71 km s" 1 Mpc" 1 ), consistent with 
WMAP measurements of the CMB power spectrum on large 
scales (Spergel et al. 2003), and quote all lengths in comoving 
units. 

2. MASS FUNCTIONS IN CDM AND IN WDM 

In this section, we briefly review the DM halo mass func- 
tions obtained in the Press-Schechter (1974; hereafter PS) and 
extended Press-Schechter (see Lacey & Cole 1993; hereafter 
EPS) formalisms, together with the modifications required to 
model structure formation in CDM models with reduced small- 
scale power. Our treatment closely follows that of BHO (which 
the reader is encouraged to consult for more details). In § 2.1, 
we describe the standard mass function, and in § 2.2, we mo- 
tivate the parameterization of a power spectrum cut-off in the 
WDM model. In § 2.3, we describe our Monte Carlo simu- 
lations needed to incorporate the additional effective pressure 
of the WDM particles with non-negligible velocity dispersion. 
Readers not interested in the modeling details of the DM halos 
can skip directly to the § 3, which describes how we associate 
GRBs with these halos. 

2.1. Halo Formation in CDM 

The mass function of DM halos in CDM models can be de- 
rived in closed form in the PS formalism. Although the PS mass 
function is in fair agreement with simulations, especially for the 
"typical" halos, it underestimates the number of rare, massive 
halos that are most relevant for our purposes; it also overes- 
timates the number of low-mass halos, when compared with 
large numerical simulations (e.g. Jenkins et al. 2001). Here 
we adopt a modified expression obtained by Sheth & Tormen 
(1999), which fits the simulated mass function to an accuracy 
of- 10%, 
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Here, dn / dM is the comoving number density of halos per unit 
mass, M is the total mass of the halo, (p) = ^MPcrit is the mean 
background matter density, 
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is the r.m.s. fluctuation in the mass enclosed within a region 
described by a top-hat filter in real space W(M) (and its Fourier 
transform Wt), and v = ^/aS c (z) / a(M), where a, p, and A are 
fitting parameters. Sheth, Mo & Tormen (2001) derive this 
form of the mass function by including shear and ellipticity in 
modeling non-linear collapse, effectively changing the scale- 
free critical over-density 5 c (z) « 1.68, obtained in the case of 
spherically symmetric collapse (Peebles 1980), into a function 
of filter scale, 
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Here b and c are additional fitting parameters (a is the same as 
in eq. 1). We use this correction to obtain the critical threshold 
S c (M,z) from 5 c (z) (in the WDM case, 5 c (z) itself is modified 
as described in § 2.3 below). For the constants in equations 
(1) and (3), we adopt the recent values obtained by Jenkins et 
al. (2001), who studied a large range in redshift and mass: a = 
0.73, A = 0.353, p = 0.175, b = 0.34, c = 0.81. 

2.2. Power Spectrum Cut-Offs 

Structure formation in WDM matter models differs from CDM 
in two main ways: (1) the free-steaming velocities of the par- 
ticles wash-out small-scale overdensities; (2) the residual par- 
ticle velocities, although they redshift away as (1+z), create 
an effective "pressure" which slows the early growth of per- 
turbations. Both effects suppress structure formation on small 
scales, but for the sake of generality we discuss the two effects 
separately. Free streaming is easily included computationally 
as a suppression of the power spectrum of fluctuations, and it 
is qualitatively similar, for example, to changes in the inflation- 
ary potential, which determines the power spectrum. Likewise, 
a red tilt of the power spectrum, n s = dlnP(k)/dlnk < 1, or 
a running of the spectral index, a = dn s jd\\\k < 0, is easily 
included in the analysis by simply modifying the power spec- 
trum in equation 2. On the other hand, the effective pressure in 
WDM models is more difficult to include computationally, and 
is specific to the WDM model; it is discussed separately in the 
next subsection. 

Free streaming manifests itself as a cut-off in the power- 
spectrum, which "flattens" a(M) for small M. This effect be- 
comes more severe as the WDM particle mass is lowered, as 
demonstrated in the bottom panel of Figure 1 for m x = 2 and 1 
keV. To produce a given energy density contribution, Q x (where 
we take tt x = 17m _ ^b), the required WDM particle mass is de- 
termined by m x n x oc Q x h 2 , where the present number density, 
n x , follows from the particle distribution function. This can be 
used (see BHO) to relate the particle mass and the r.m.s. veloc- 
ity dispersion, v„ 
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where g x is the effective number of degrees of freedom of WDM. 
The usual assumption of a fermionic spin-^ particle yields g x = 
1.5. This modifies the CDM power spectrum (which we take 
from Eisenstein & Hu 1999) by multiplying it with the square 
of a transfer function (Bode, Ostriker, & Turok 2001): 

T x {k) = {\ + {ekR c ) 2v y^ v , (5) 
with parameters e = 0.361, r) = 5, v = 1.2. The power spectrum 
is reduced to half its value in CDM at k = 1/R C , where the cut- 
off scale, R c , is given by 

'O h 2 \°' 15 / e \-0.29/ m N -1.15 



where R c is in comoving Mpc. Here, we only consider parti- 
cles with fermionic spins, i.e. g x = 1.5, but all results can be 
scaled for arbitrary values of g x using the equation above. The 
corresponding mass scale, M c = (A/3)ttRI(p), is 
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We will use R c and m x interchangeably, with equation (6) relat 
ing them. 
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FIG. 1 . — Top panel: critical overdensity threshold as a function of scale, 
evaluated for z = 10. The dotted curve represents the scale-free critical over- 
density 8 c (z). The solid curve corresponds to CDM, the long-dashed curve 
correspond to a WDM particle mass of m x = 2 keV , and the short-dashed 
curves correspond to m x = 1 keV . All curves incorporate the correction from 
equation (3), except for the dotted curve. Bottom panel: the r.m.s. mass fluc- 
tuation in top-hat filters of mass-scale M. As above, the solid curve corre- 
sponds to CDM, the long-dashed curve corresponds to m x = 2 keV (or power- 
spectrum cut-off scale R c = 0.087 Mpc), and the short-dashed curve corre- 
sponds tom x -l keV (or R c = 0.193 Mpc). 



2.3. Effective Pressure of WDM Particles 

As mentioned above, structure formation in WDM models 
is further suppressed by the residual velocity dispersion of the 
WDM particles, which delay the growth of perturbations. We 
use the results of BHO, who made an analogy with an ideal gas 
whose temperature corresponds to the velocity dispersion of the 
WDM. In this case, the pressure delays the collapse, and can be 
effectively included in the mass function computed in the EPS 
analysis, by raising the critical linear extrapolated overdensity 
threshold at collapse 5 c (z). BHO computed the critical overden- 
sity by following the collapse of spherical perturbations, us- 
ing a one-dimensional, spherically symmetric Lagrangian hy- 
drodynamics code, originally developed by Thoul & Weinberg 
(1995). The top panel in Figure 1 shows 5 C (adopted from BHO) 
as a function of M, at the fixed redshift z = 10. The dotted curve 
represents the scale-free critical overdensity 6 c (z) arising from 
spherical collapse. The solid curve corresponds to CDM, the 
long-dashed curve corresponds to WDM with m x = 2 keV , and 
the short-dashed curve corresponds to m x = 1 keV. The thresh- 



4 



old is found to increase sharply below the mass scale 

which can be shown (see BHO) to correspond to an analog of 
a "Jeans mass", i.e. the scale of the objects whose collapse is 
significantly delayed by the pressure. 
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FIG. 2. — The fraction of the total mass collapsed into halos of mass M 
or higher, F(> M,z), as a function of redshift. The shaded region shows the 
collapsed fraction in CDM, with a range of low-mass cutoffs corresponding 
to virial temperatures 300 K < T ylI < 10 4 K. The other curves correspond to 
WDM particle masses of m x = 3.0 keV, 2.5 keV, 2.0 keV, 1.5 keV, 1.0 keV, 
0.5 keV, from top to bottom, and do not assume any low-mass threshold (ef- 
fectively, M = 0). The cutoffs used in the CDM case would leave the WDM 
results essentially unchanged (see discussion in text). 

Since S c (M,z) is function of scale in WDM models, one can 
not obtain WDM mass functions from the standard EPS anal- 
ysis, which makes use of the symmetry in the random walk 
trajectories of 5 vs. M about a fixed threshold (Lacey & Cole 
1993). Instead, we compute the mass functions using Monte- 
Carlo simulations. We generate random realizations of trajecto- 
ries 5(M), as the scale is decreased fromM <~ oo, and generate 
the histograms of the scales at which the trajectory first crosses 
the 6 C (M, z) threshold. Each step in the random walk, AS(M), is 
Gaussian distributed with a variance of <7 2 te _(M, AM) = cr 2 (M) - 

a 2 (M- AM). When constructing such a random walk, one must 
be careful to use steps small enough such that as the smoothing 
scale M is decreased by AM, the likelihood that the 5 c (M,z) 
threshold is crossed between M and M - AM is small. The 
physical reason for this is the so-called "cloud-in-cloud prob- 
lem": to ensure that we do not step over a collapsed halo as we 
decrease the smoothing scale M (i.e. that our 8(M) trajectory 
has not gone above S c (M,z) and then dropped below it again 
within the range AM). We use an adaptive step-size, set so 
that the barrier 5 c (M,z) is at least 7cr step away from S(M), with a 
minimum resolution of AM = M/100. Formally, defining A 7 M 



such that (6 c (M,z)-S(M)) I er ste p(M , A 7 M) = 7, our step size is 
AM = MAX[M/100, A 7 M]. We find that these parameters effi- 
ciently reproduce the standard EPS mass function in the CDM 
case to an accuracy of a few percent. 

We also find that a starting mass for the random walk tra- 
jectories as small as <~ 10 16 M Q is sufficient to obtain accurate 
mass functions at z = 0. The starting mass can be decreased 
as redshift increases, since the characteristic mass that is col- 
lapsing gets smaller as redshift increases. A starting mass of 
— 10 12 M Q is sufficient at z = 15. 

The number of simulated trajectories required to obtain ac- 
curate mass functions is a strong function of redshift. This is 
to be expected, since virialized structure becomes very rare at 
high redshift. We show this effect in Figure 2, where we plot 
the fraction of the total mass collapsed into halo of mass M or 
higher, F(> M,z), as a function of redshift. The shaded region 
shows the collapsed fraction in CDM, with a range of low-mass 
cutoffs corresponding to virial temperatures 300 K < T wv: < 10 4 
K (see discussion in § 3 below). The other curves correspond 
to WDM particle masses of m x = 3.0 keV, 2.5 keV, 2.0 keV, 1.5 
keV, 1.0 keV, 0.5 keV, from top to bottom, and do not assume 
any low-mass threshold (effectively, M = 0). Introducing the 
same two low-mass cutoffs as in the CDM case would leave 
the WDM results essentially unchanged, since the power is al- 
ready strongly suppressed in excess of these cutoffs at the high 
redshifts where these low mass halos would dominate the col- 
lapse fraction (see discussion below, or Fig. 7 in BHO). As can 
be seen, simulating high redshift mass functions for small parti- 
cle masses can be prohibitively computationally expensive. For 
example, we find that accurate mass functions for m x < 1 keV 
at z > 15, require <; 10 9 Monte-Carlo runs, as less than one in 
a million S(M) trajectories crosses S c (M,z). 
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FIG. 3. — Cumulative mass functions (number of halos with masses greater 
than M per Gpc 3 ) at redshifts z=10 (top panel) and z=15 (bottom panel). The 
solid curves show mass functions in CDM models; the dotted curves corre- 
spond to WDM models with m x = 2 keV and 1 keV (from top to bottom); the 
dashed curves are mass functions in the same two WDM models, but incorpo- 
rating only a power spectrum cut-off (§ 2.2), and no WDM pressure (§ 2.3). 
The WDM cumulative mass functions were created with 10 9 Monte-Carlo re- 
alizations of 5(M) trajectories. 
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Figure 3 shows sample cumulative mass functions at red- 
shifts z=10 {top panel) and z=15 {bottom panel). The solid 
curves correspond to CDM; the dotted curves correspond to 
WDM models with m x = 2 keV and 1 keV (from top to bottom). 
The dashed curves are mass functions for the same two WDM 
models, but incorporate only a power spectrum cut-off, ignor- 
ing the effective pressure of WDM. The WDM mass functions 
were created with 10 9 Monte-Carlo runs as explained above. 
(Note that our results are slightly different from those of BHO; 
this is due to a small correction to BHO's derivation of R c and 
the corresponding <r{M).) 

As can be seen from Figure 3, including the pressure term in 
the WDM models further suppresses the number of halos rela- 
tive to the models which include only the power spectrum cut- 
off. Furthermore, the relative importance of the pressure term 
increases with increasing redshift. As anticipated, the overall 
differences between the WDM and CDM mass functions in- 
crease toward higher redshift. This is because in the early uni- 
verse, the characteristic scale of collapsing and virializing halos 
was smaller, and closer to the cut-off scales discussed above. 
These large differences will aid in discriminating between mod- 
els with different power spectrum cut-off scales. 

3. THE EVOLUTION OF THE GRB RATE WITH REDSHIFT 

In this section, we describe a model for the expected evolu- 
tion in the rate of all GRBs, as well as the fraction that can be 
detected by Swift. As is commonly done in the context of pre- 
dicting the rate of GRBs at high redshift (Bromm & Loeb 2002; 
Choudhury & Srianand 2002), we assume that the GRB rate 
density (the number of GRBs per unit time per unit comoving 
volume), pgrb{z), is proportional to the global star-formation 
rate density, 

Pgrb(z) ~ K p\{z) , (9) 

where K is the proportionality constant in units of Mq, and 
p\{z) is the stellar mass produced on average per unit comoving 
volume per unit time. First we discuss the evolution of the star 
formation rate, and then the normalization of the correspond- 
ing GRB rate. In § 3.1.2 we discuss the uncertainties of our 
approach. 

3.1. Star-Formation Rate in Halos 

We estimate the global star-formation rate (SFR) density at 
redshift z as 

J Mmm J oo dMdz' 
where e* is the efficiency parameter for the conversion of gas 
into stars, M dM dn{> M,z')/dM is the mass density con- 
tributed by halos with masses between M and M + dM at red- 
shift z, t{z) is the age of the universe at redshift z, and P{t) is the 
probability per unit time that new stars form in a mass element 
of age r = t{z)-t{z'). We adopt the fiducial value of e» =0.1 
(see, e.g., Cen [2003]), but note that our results are insensitive 
to this value since we normalize the coefficient K in order to 
match our total GRB rate with observations (the important as- 
sumption is only that e* is constant; this assumption is probably 
conservative, as discussed below). We also note that the simple 
time-derivative of the halo mass-function in equation (10), in 
general, contains a contribution from mergers between halos, in 
addition to the formation of new halos; see, e.g., Sasaki 1994 
for a discussion. However, the integral in equation (10) is sen- 
sitive only to the total virialized mass above mass M m j n , so that 



this ambiguity should not affect our results (although we note 
that the formation of new halos dominates in our case since the 
relevant halos are above the non-linear mass-scale). 

The minimum mass, M m ; n , depends on the efficiency of gas 
cooling and collapsing into a dark matter halo. Without molec- 
ular hydrogen, M m [ n corresponds to a halo with virial temper- 
ature, T vir - 10 4 K; with H 2 , T vir - 300 K (Haiman, Abel, & 
Rees 2000; we use the conversion between halo mass and virial 
temperature as given in Barkana & Loeb 2001). The amount 
of H2 present in the early universe is uncertain, so below we 
present results for both r v ; r > 10 4 K and T vn > 300 K. We note 
that only the CDM mass functions are sensitive to these cooling 
thresholds (see Fig. 4). In our models which incorporate small- 
scale power suppression, as the characteristic collapse scale 
approaches the cooling cut-off, T wu , the cooling cut-off has 
already become smaller than the power spectrum cut-off, R c . 
Hence, we do not distinguish between models with r v ; r > 10 4 
K and T vn > 300 K for the WDM particle masses presented be- 
low. 

In the presence of an ionizing background, the cosmological 
Jeans mass, above which gas can collect in DM halos, is in- 
creased, corresponding to virial temperatures 10 4 K (Rees 
1986; Efstathiou 1992). Early work on this subject (Thoul & 
Weinberg 1996) suggested that an ionizing background would 
completely suppress star formation in halos with circular ve- 
locities v c ; rc <, 35 km/s, and partially suppress star-formation 
in halos with 35 km/s <, v C i rc <^ 93 km/s. More recently, Di- 
jkstra et al. (2004) found that such a suppression is likely to 
be countered by strong self-shielding for z ^ 3. For con- 
creteness, we completely suppress star-formation in halos with 
T v i r < 10 4 K for z < 7 (assumed to correspond to the redshift 
of reionization, when the background radiation is established; 
e.g. Mesinger & Haiman 2004), and increase the cut-off to 
Vcirc < 55 km/s (r vir < 1.1 x 10 s K) for z < 3 (when the inter- 
mediate halos with 10 km/s < v c i rc < 55 km/s are assumed to 
cease self-shielding). As can be seen in Figure 4 below, these 
assumptions only weakly affect our results, because most of 
the contribution to p\{z) at these redshifts comes from more 
massive halos. Furthermore, we already bracket a large range 
of expected GRB distributions by presenting results for both 
T vir > 10 4 K and T vil > 300 K. 

Finally, we describe two different functional forms for the 
stellar formation probability density P{t). 

3.1.1. Instantaneous Star-Formation 

The simplest assumption is to adopt P{t) to be a Dirac delta 
function, 

P{t) = S{t). (11) 

This essentially assumes that a fraction e*(£!/,/f2 M ) of the cur- 
rently virializing mass is instantaneously converted into stars, 
and that previously virialized mass does not contribute to p*(z). 
With this assumption, equation (10) becomes 

fl b f°° i/r dz d 2 n{>M,z) ... .... 

/9*(z) = e»— / M— — — — dM . (12) 

J Mmin dt dzdM 

Since the WDM mass functions are obtained with computa- 
tionally expensive Monte-Carlo simulations, we only present 
their results for this simple model; results for mass functions 
incorporating a power spectrum cut-off are presented for both 
P{t) models. In numerically calculating derivatives of the mass 
functions in WDM (An(> M,z)/Az), we use z step-sizes of 
Az ~ 0.2 - 0.5 (at high redshifts An(> M,z)/Az is a flatter 
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function of z, so good accuracy can be achieved even with a 
larger Az step-size). The resulting SFR densities are shown 
in the top left panel of Figure 4. The solid curves assume 
r v i r > 10 4 K, and include only the effect of the power-spectrum 
cutoff, assuming cut-off scales of R c = Mpc (m x — > oo; i.e. 
standard CDM), R c = 0.193 Mpc (m x = 1 keV), and R c = 0.087 
Mpc (m x = 2 keV), from top to bottom. The dashed curves as- 
sume r v ; r > 300 K, with the same cut-off scales. The dotted 
lines include the additional effect of the effective pressure of 
WDM particles with a mass of m x = 2 keV and m x = 1 keV, from 
top to bottom. Since power is strongly suppressed on scales in 
excess of r vir > 10 4 K for these WDM particle masses cut-off 
scales, results for T vu > 10 4 K and r v ; r > 300 K are virtually 
the same. In CDM models, the range of virial temperatures 
result in a wider range of expected SFR densities, highlighted 
by the shaded area. Our results are within the uncertainties of 
existing observational estimates of the SFR density in the cur- 
rently available, low-redshift (z <^ 2) regime (e.g. Schiminovich 
et al. 2004). Our SFR density increases out to z <; 6 (and to 
still higher redshift in the T vn > 300 K case), and is consistent 
with recent estimates at 2 <, z < 6 from the Hubble Ultra Deep 
Field (Bunker et al. 2004), the GOODs surveys (Giavalisco et 
al. 2004) and the FORS Deep Field on the VLT (Gabasch et 
al. 2004), after they incorporate a factor of 5-10 increase in the 
SFR (Adelberger & Steidel 2000) due to dust obscuration. Note 
that the peak of the SFR density falls between 3 ^ z ^ 12 in our 
range of models, and that the SFR remains significant (exceed- 
ing its present day value) out to redshifts as high as z <; 25. 

We stress that instant star-formation is the most conserva- 
tive model in setting observational limits to a power spectrum 
cut-off, as it predicts the largest number of high-redshift GRB 
detections. Any other, more realistic form for the stellar forma- 
tion probability density, P(t), will delay GRB events, thereby 
smearing the event rate toward lower redshifts, and decreasing 
the likelihood of detecting high-redshift GRBs, especially in 
the presence of small-scale power suppression (see § 4.1). 

3. 1 .2. Exponential Approximation for P(t) 

We alternatively assume that stellar formation occurs on a 
time-scale corresponding to the dynamical time, td yn ~ (Gp)~^ 2 , 
(Cen & Ostriker 1992; Gnedin 1996): 



P(t(z)-t(z'))- 



t(z)-t(z') 



-exp 



l dyn 



tjz)-t(z') 

$dyn 



(13) 



as shown in the top left panel, except the WDM models with 
pressure (dotted curves) are not included. As expected, intro- 
ducing a finite width to the stellar formation probability density 
preferentially suppresses high-redshift star-formation. 

We remark that the sharp drop at z = 7 in the T wu > 300 K, 
CDM curve is an indicator of the assumed sudden reionization 
(Bromm & Loeb 2002; Choudhury & Srianand 2002). This fea- 
ture is a direct prediction of reionization models, and is indeed 
likely to be sharp (see, for example, Cen & McDonald 2002 
for a discussion, and a claim for detecting a related feature in 
the Lyman a transmission spectrum of distant quasars). Given 
a sufficient number of GRBs, the shape and location of such a 
feature can be mapped out, and used as a probe of the epoch 
of reionization. We postpone a detailed exploration of the de- 
tectability of such a feature to a future paper. 



where p(z) « A c p crit (z) is the mean mass density interior to col- 
lapsed spherical halos (e.g. Barkana & Loeb 2001), and A c is 
obtained from the fitting formula in Bryan & Norman (1998), 
with A c = 187r 2 178 in the Einstein-de Sitter model. 

Since there is no unique, physically-motivated, and self-consistent 
way to track individual mass elements and halo mergers in the 
EPS formalism, assigning an age, r = t(z)-t(z'), to each mass 
element is somewhat arbitrary. The problem arises since two 
neighboring mass elements which are part of the same halo, 
can in EPS be flagged as belonging to two different halos with 
different masses (Somerville & Kollatt 1999; see also Sheth & 
Pitman 1997 and Benson et al. 2004). Here, we assume that 
a mass element "starts its clock" (t = 0) when it first becomes 
part of a halo with mass M > M m ; n , and that it carries around 
that clock through any future mergers without resetting it. This 
ambiguity is bypassed in § 3.1.1, since P(t) is assumed to be a 
delta function. 

We present our SFR densities for this model in the top right 
panel of Figure 4. The curves correspond to the same models 
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FIG. 4. — Top Panel: SFR densities assuming instantaneous star-formation 
in a collapsing halo (left panels) and a finite exponential spread in star- 
formation times (right panels). The solid curves assume a minimum virial 
temperature 7" v j r = 10 4 K of a halo for star-formation, and include the effect 
of power-spectrum suppression below various cut-off scales: R c = Mpc 
(m x -> oo; i.e. standard CDM), R c = 0.193 Mpc (m x = 1 keV), and R c = 
0.087 Mpc (m x = 2 keV), from top to bottom. The dashed curves assume a 
minimum temperature, r v ; r = 300 K, with the same cut-off scales. The dotted 
curves correspond to models that include the additional effect of the effective 
pressure of WDM particles with mass m x = 2 keV and m x — \ keV, from top 
to bottom. The shaded area highlights the expected range in SFR densities in 
CDM models with minimum virial temperatures 300K < r v ; r < 10 4 K. Bottom 
Panel: The fraction of all GRBs that originate at redshifts higher than z in the 
models corresponding to the top panels. 



3.2. The GRB Rate Associated with Star-Formation 

In order to constrain the proportionality constant K by match- 
ing the predicted and observed GRB rates, a luminosity func- 
tion (LF) for the GRB population needs to be assumed. Ob- 
servations have shown that GRBs are far from being standard 
candles (and even if the total GRB energy has a nearly universal 
value, the burst luminosity will vary; Frail et al. 2001). How- 
ever, the determination of their intrinsic LF has been hampered 
so far by the lack of a sufficiently large sample with detected 
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redshifts. On the other hand, fits to the observed flux distribu- 
tion suffer from a degeneracy between the LF and the SFR that 
prevents an independent determination of both quantities. 

On the theoretical side, simulations of jets in the collapsar 
model (MacFadyen & Woosley 1999) have shown that the jet 
energy (and hence luminosity) is a decreasing function of the 
viewing angle 9 relative to the jet axis, triggering studies of 
structured jets (Rossi et al. 2002; Zhang & Meszaros 2002). 
Here, following Rossi et al. (2002), we assume that the lumi- 
nosity has a form L = L m \ n (9j et /9) 2 for 9 COK < 6 < f?j et , where 
9 cole defines the core of the jet (within which the luminosity is 
constant), and #j et defines the outer edge of the jet (above which 
the luminosity drops to zero). Current observations have shown 
that 6» core <, 0.06 rad, and jet ^ 0.6 rad (Bloom et al. 2003). 
Our results are insensitive to the precise value of 6> core below 
the observed minimum value, and therefore we set it to zero for 
simplicity. The GRB rate at z = 0, is, however dependent on 
#j et (or equivalently, the inferred minimum luminosity L m i n ; see 
also Guetta et al. 2004). 




FIG. 5. — The fraction of all GRBs that are brighter in our model luminos- 
ity function than the detection thresholds of BATSE (dashed curves) and Swift 
(solid curves). The lower/upper curve in both cases corresponds to jet angles 
of Sjei = t/2 and 0.6 rad, respectively. In each case, we determine the mini- 
mum burst luminosity L m j n so that the flux distribution predicted in our CDM 
model with a threshold virial temperature T vir = 10 4 K best fits the distribution 
observed by BATSE. 

For each of the models for the SFR described in the previ- 
ous section, we determine the constant L m i n by finding the best 
fit between the theoretical and the observed flux distribution of 
bursts (unnormalized; i.e. using only the shapes of the distri- 
butions). In this fitting procedure, we fix the jet angle to be 
f5j et = 7r/2 (but consider the alternative #j et = 0.6 rad, the low- 
est value directly inferred from observations so far; see below). 
In practice, a distribution of viewing angles is inferred for the 
theoretical GRBs (assuming random jet orientations), and com- 
pared to the distribution dn / d9 inferred from the bursts detected 
by BATSE (see Perna et al. 2003 for details of the analysis). 
The proportionality constant K is then separately determined 



by imposing that the total number of bursts J d6(dn/d6) 
above the BATSE sensitivity (~ 0.25 ph sec" 1 cm" 2 ) be 667 
yr" 1 . Once these two model parameters are determined, the 
number of bursts observed at redshifts greater than z, over a 
time interval, Af b s , and solid angle, Afi, 



An . 

N(> z) = —At obs 
4tt 



dz 



, Pgrb(z') dV(z') 
(1+zO dz' Jo 



x(z') 



P(8)d6, 
(14) 



can be found for the Swift sensitivity (F\i m <~ 0.04 ph sec" cm ). 
We adopt a solid angle Ail = 0.34 sr, for the fully coded detec- 
tor area of Swift. In the above equation P(9)d9 = sin 9d9 is the 
probability of viewing a randomly oriented GRB at an angle 
from the jet axis between 9 and 9 + d9, the factor 1 /( 1 + z) ac- 
counts for time dilation, and dV(z)/dz is the comoving volume 
in our past light cone per unit redshift, 



dV(z) 



■■Aire 



4(z) 



(15) 



dz l+z 

where c is the speed of light, and <fz.(z) is the luminosity dis- 
tance. The maximum angle 9 m . dx at which a burst can be de- 
tected at a given redshift is found by numerically inverting the 
equation 

F — ^(^max) Hft\ 

Hm ~ 4ndt(z)(l + zr' i b) 

where a is the power-law spectral index of the photon number 
spectrum of the burst (see Lamb & Reichart 2000 for details). 

For illustration, in Figure 5, we show the fraction of all bursts 
at redshift z that would be detectable by BATSE (dashed curves) 
and by Swift (solid curves). The lower/upper curves correspond 
to jet angles of 9 jsl = ir/2 and 0.6 rad, respectively. In our fidu- 
cial model, Swift detects all GRBs out to z ~ 1, and <~ 1% of 
GRBs atz~ 10. 

The value of the normalization constant, K, inferred from the 
BATSE data, is dependent on the value of L m j n . For example, in 
the CDM models with instantaneous star-formation and r v j r > 
10 4 K, for our fiducial choice of #j et = ir/2, we find the best- 
fitting value of L m ; n = 1 .5 x 10 55 ph/sec, and the local GRB rate 
of Pgrb(z = 0) = 6 Gpc" 3 yr" 1 , while for 6>j et = 0.6 rad, we find 
instead L m ; n = 10 56 ph/sec, and pgrb(z = 0) = 1 Gpc" 3 yr" 1 . 

Nevertheless, we find that our subsequent predictions for the 
number of bursts at high redshift by Swift are only mildly de- 
pendent on the choice of 6>j et . This is because as L m ; n is in- 
creased, the normalization factor K decreases. This tends to 
compensate for the increase in the observable number of bursts 
due to the fact that there is a larger fraction of bursts at higher 
luminosities. As can be seen in Figure 5, at high redshift only 
the bright power-law tail of GRBs are detectable, and the two 
effects nearly cancel. We note that the slope of N(> z) depends 
on the power-law slope of the LF. Our model has dn/dLct L" 7 
with 7 = 2 in the high luminosity tail - for comparison, by fit- 
ting a single power-law to dn jdh, Firmani et al. (2004) find a 
shallower slope with 7 <~ 1.6. If this slope is accurate for high 
luminosities, it would somewhat increase the number of high- 
Z GRBs. Of course, once Swift data has been gathered, it will 
provide new and independent constraints on the GRB luminos- 
ity function, facilitating more accurate estimates of the high-z 
GRB rate. 

3.3. Uncertainties in GRB Redshift Evolution 

Our most important model assumptions are that e* and K are 
constant. Observations of nearby dwarf galaxies (Taylor et al. 
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1999; Walter et al. 2001) yield a range of efficiencies w 0.02 
- 0.08. These values correspond to the low halo mass scales 
which form at high redshifts, although there is no clear direct 
correspondence between the first generation halos and these lo- 
cal dwarfs. Numerical simulations of metal-free star-formation 
at high redshift (Abel, Bryan, & Norman 2002; Bromm, Coppi 
& Larson 2002) suggest that the first generation of stars form in 
minihalos with even lower efficiencies, <^ 0.01. Such a re- 
duction in minihalo star-formation (r v ; r < 10 4 K) suggests that 
at high redshifts, the true GRB rate is closer to our T v i r > 10 4 K 
curve. 

On the other hand, since GRB progenitors involve high-mass 
stars (MacFadyen & Woosley 1999), and the stellar initial mass 
function (IMF) is expected to be more top-heavy in the early 
universe, K might increase at high redshifts, increasing pgrb(z)- 
We note that the minimum mass of stars that loose their hydro- 
gen envelope and may thus produce GRBs can also strongly de- 
pend on metalicity. Heger et al. (2003) find that this minimum 
mass increases from w 30M Q at solar metalicity to w 100M Q 
at zero metalicity (see their fig. 3). However, the trend is still 
more likely to be for K to increase toward high redshift, be- 
cause > 30M Q stars are relatively rare in local galaxies with a 
Salpeter IMF, while the simulations of Abel et al. (2002) and 
Bromm et al. (2002) suggest that metal-free star formation at 
high redshift may produce exclusively > 100M Q stars. Thus, 
since they are expected to act in opposite directions, the red- 
shift evolutions of and K might compensate somewhat for 
each other. The rather high angular momentum of the collaps- 
ing star required to produce a GRB is more easily achieved if 
the star is in a binary, and binaries are found to be more fre- 
quent at low metallicity (i.e. high z; Fryer et al. 1999). It is, 
however unclear if binaries do form in the first generation of 
truly metal-free halos. Simulations following direct cosmolog- 
ical initial conditions by Abel, Bryan & Norman (2002) find 
a single star with no further fragmentation; however, different 
simulations with somewhat more artificial initial conditions (a 
rotating cylinder), find efficient binary-formation (Saigo et al. 
2004). 

It is conceivable that at high redshifts, a WDM universe could 
mimic the GRB rates of a CDM universe, by compensating for a 
loss of small-scale power with higher efficiencies of GRB pro- 
duction (e.g. a higher K value at high redshifts), and visa-versa. 
However, we note that such a redshift evolution is unlikely to 
be sharp enough to significantly impact our conclusions. Even 
a power-law evolution of K would be insufficient to compen- 
sate for the exponential suppression of small-scale power in our 
models. For example, from Figure 6 we note that in order for 
GRB rates in WDM models with m x w 2 keV to match GRB 
rates in CDM models at z > 10, the product e*K would need 
to be a factor of 10 higher than the overall average value deter- 
mined from the lower-redshift BATSE sample. For m x w 1 keV, 
e*K would have to be <~ 100 times larger at z > 10. These dif- 
ferences only increase with increasing redshift, and if such an 
increase is present, it can, in principle, be detected by studying 
the shape of the GRB redshift distribution function (see § 4.2). 

Recently, Wise & Abel (2004) calculated primordial super- 
novae (SNe) rates with a semi-analytic analysis of feedback 
mechanisms and evolution of primordial stellar environments, 
constraining their results with the measured WMAP optical depth 
to electron scattering, r e = 0.17. With ab-initio knowledge of 
GRB progenitors, an analogous analysis could be performed to 
predict the GRB rate as a function of redshift, replacing the pro- 



portionality constant, K, assumed in this paper. In practice, our 
current knowledge of the physics of GRB progenitors and of the 
various feedback processes in the early universe is highly un- 
certain, and such an approach would introduce additional free 
parameters. Nevertheless, this approach could be useful in the 
future, since some of the free parameters may be independently 
constrained. For example, in addition to GRBs, early star- 
formation will be accompanied by observable SN explosions 
(e.g. Miralda-Escude & Rees 1997), the production of heavy 
elements (e.g. Haiman & Loeb 1997), remnant stellar black 
holes (e.g. Volonteri et al. 2003; Volonteri & Perna 2005), 
and the reionization of the IGM. Observations of these effects 
can further reduce the uncertainties associated with early star- 
formation history and GRB progenitors, and thus ultimately re- 
duce uncertainties on the expected primordial GRB event rates 
(for a discussion of such future constraints, see Wise & Abel 
2004, who also emphasize that SNe observations can constrain 
many properties, such as mass, luminosity, metallicity, and red- 
shift, of the progenitors). Finally, early "mini-galaxies" may 
also be directly detectable (Haiman & Loeb 1997) by the James 
Webb Space Telescope (JWST). 2 

4. CONSTRAINTS ON STRUCTURE FORMATION MODELS 

The model outlined in the previous sections can be used to 
compute the evolution of the GRB rate with redshift, as well 
as the flux distribution of the bursts, allowing us to incorporate 
the detection threshold of Swift. In this section, we present con- 
straints on WDM models first from the total number of GRBs, 
and then a potentially less model-dependent constraint from the 
distribution of a luminosity-limited sub-sample of the bursts. 
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FIG. 6. — Expected Swift detection rates of GRBs occurring at redshifts 
greater than z. The curves correspond to the same models as shown in Figure 4, 
and the horizontal line denotes the detection threshold of 1 GRB/yr. The top 
panel assumes instantaneous star-formation; the bottom panel assumes a finite 
exponential spread in star-formation times. 



4. 1 . Absolute GRB Detection Rates 
2 See www.jwst.nasa.gov 
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The most straightforward constraints follow from the total 
Swift burst detection rates. In Figure 6, we show the cumulative 
Swift GRB detection rates produced by our models. The curves 
correspond to the same models as shown in Figure 4 above, 
with the shaded area enclosing the range of expected distribu- 
tions in CDM models for different limiting virial temperatures 
for star-formation in early halos. The top panel assumes instan- 
taneous star-formation, and the bottom panel assumes a spread 
of star-formation times as discussed above. The horizontal line 
denotes a detection rate of 1 GRB per year, roughly the lower 
limit for detections in Swiff s two-year nominal operation. 

Figure 6 shows that detections of bursts originating at red- 
shifts larger than z ~ 15 are improbable even in CDM models. 
Constraints on WDM models can be simply read off this figure 
by noting the redshifts at which the curves corresponding to a 
given model intersect the horizontal line. For example, for a 
WDM particle mass of m x = 1 keV, a detection rate of one z > 7 
GRB per year is expected. That number drops to one z > 8 
GRB per 10 years, which is unlikely to be detectable. A de- 
tection of a single z k 10 GRB by Swift would rule out models 
incorporating WDM or equivalent power spectrum cut-offs of 
m x <,2 keV(R c k 0.087 Mpc), and would therefore significantly 
improve the existing constraints. 

These constraints on models with an exponential cut-off in 
the small-scale power spectrum are summarized further in Ta- 
ble 1 . The low end of the quoted redshift range in each case in 
the top five rows corresponds to the exponential star-formation 
and the high end to instantaneous star-formation (the models 
that include the effective pressure of WDM particles, listed in 
the bottom three rows, show only the latter). The right-most 
column in the table shows the redshifts beyond which 10 GRBs 
per year would be detectable. 

The constraints shown in Table 1 have been derived assum- 
ing a fixed cosmological model, varying only the cut-off scale 
for the power spectrum. However, in principle, the effects of 
the cut-off can be mitigated by changing other parameters. In 
particular, the GRB rate depends exponentially on the normal- 
ization of the power-spectrum a%. In order to quantify the 
robustness of our results, we have re-computed the redshifts 
corresponding to 1 and 10 GRBs/yr in the CDM model with 
r v ; r > 10 4 K and in the model with the power-spectrum cutoff 
at R c = 0. 193 Mpc, but using a high value of erg = 1 rather than 
our fiducial choice of a& = 0.85. We have found that this in- 
creases the redshifts by Az < 2; in particular, the R c = 0.193 
Mpc model could still be ruled out by observing 1 GRB/yr at 
zk9.1. 



Table 1 

The redshift, in several WDM models, beyond which 
the Swift satellite would detect 1 or 10 GRBs/year. 



The corresponding constraints in models with a red tilt or a 
running of the scalar index are summarized in Table 2. Overall, 
in the CDM model with the more stringent limit of r v i r > 10 4 K 
for star-formation, the effect of introducing a tilt with n s = 0.9, 
or a running of the index with a = dn s /dlnlc = -0.05 is compa- 
rable to that of the exponential cut-off in the power spectrum 
due to a m x ss 2 keV WDM particle: the GRB rate falls below 
1/yr at z k 8. However, as the table shows, these constraints are 
less robust than those for WDM. The suppression of the power 
spectrum is a shallow function of scale, and, unlike in the case 
of the exponential suppression, the GRB rate can be increased 
by either a relatively modest change in the normalization erg, or 
by allowing star-formation in smaller halos down to r v ; r = 300 
K. The table suggests that a detection rate of 1 GRB/yr even 
at z k 14 cannot be used to distinguish a tilt of n s w 0.9 from 
models with n s = 1 but with a lower r v ; r > 300 K and a higher 
as. However, a detection rate of 1 GRB/yr at z k 10 would lead 
to a relatively robust upper limit on the running of the spectral 
index, a k -0.05. 



Table 2 

The redshift in models with a red-tilted or running 
scalar index, fixed Swift GRB rates, as in Table 1. 



Model 



1 GRB/yr 



10 GRBs/yr 



Red Tilt (n s = 0.9) 



r vir > io 4 k 


zk 7.7- 


9.3 


zk 4.1 


-5.0 


T vir > 300K 


zk 10.5 


-12.6 


zk 4.8 


-6.3 


a s = l (r vir > io 4 K) 


zk 9.1- 


11.1 


zk 4.8 


-5.8 


cr 8 = 1 (r vir > 300K) 


zk 12.2 


- 14.7 


zk 5.8 


-7.3 



Running (a = -0.05) 



r vir > 10 4 K 
r vir > 300K 

a s = i (r vir > io 4 K) 

o- 8 = l (r vir > 300K) 



zk 7.4-8.9 

zk 8.8-10.5 

zk 8.7-10.5 

zk 10.2-12.3 



zk 4.0-4.9 

zk 4.3-5.4 

zk 4.6-5.6 

zk 5.1-6.4 



4.2. Redshift Distribution of a Luminosity-Limited Sample 

A method that is less dependent on models of the GRB LF, 
and that could, in principle, discriminate against WDM models 
and models with a power spectrum cut-off is to construct a red- 
shift distribution of a luminosity-limited sub-sample of the ob- 
served bursts. Under the relatively weak assumption that the LF 
of GRBs does not evolve with redshift (but without any other 
assumptions about the LF), this redshift distribution would be 
proportional to the distribution of all bursts. In the bottom pan- 
els of Figure (4), we plot the fractional distribution of all bursts, 



Model 


1 GRB/yr 




10 GRBs/yr 


CDM (r vir > 10 4 K) 


zk 9.40- 


-11.6 


zk 4.8 


-5.9 


CDM (r vir > 300K) 


zk 13.4- 


- 16.0 


zk 6.4 


-7.6 


P(k) suppr (R c £ 0.193Mpc) 


zk 6.3- 


7.5 


zk 3.5 


-4.3 


P(k) suppr (R c <> 0.087Mpc) 


zk 7.8- 


9.3 


zk 4.2 


-5.1 


P(k) suppr (R c <, 0.055Mpc) 


zk 8.4- 


10 


zk 4.5 


-5.4 


WDM(m^ IkeV) 


zk 6.7 




zk 3.9 




WDM (m x k 2keV) 


zk 8.5 




zk 4.8 




WDM(m^ 3keV) 


zk 9.5 




zk 5.2 





5 , , Pgrb(z') dV(z') 
dZ (1+zO dz' X 



f 

Jo 



dz 



,Pgrb(z') dV{z') 
(1+zO dz' 



(17) 

This quantity is then independent of the value of the normaliza- 
tion K and of the shape of the GRB LF discussed above. 

Such intrinsic event rate distributions can be constructed by 
tracking the redshift distribution only for the subset of GRBs 
above a minimum luminosity that allows one to see them up to 
a given high redshift, and neglecting GRBs which are fainter 
than this threshold. By construction, all of the remaining GRBs 
are luminous enough to be detectable up to the chosen redshift 
and so trace out the intrinsic event rates, / (i.e. independently of 
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the shape of the LF and K). The high sensitivity of Swift makes 
this still a sizable sample. For example, in the case of the stan- 
dard CDM model, we find that w 30% of the bursts would still 
remain by choosing only the ones above the minimum luminos- 
ity that allows for their detection up to z ~ 10. Likewise, w 20% 
of bursts would constitute the sample with the minimum lumi- 
nosity that allows for their detection up to z ~ 15. The bottom 
panels of Figure 4 reveal that <~ 50% of these remaining bright 
GRBs would be at redshifts z^5, where the curves start show- 
ing significant deviations from one another. In the CDM case, 
this would translate to a sample of ~ 20 high-redshift (z 5) 
GRBs with which to trace out the redshift distribution in the 
bottom panels of Figure 4. 

Finally, we note that even for this high tail of the LF, the 
efficiency with which ultra-high redshift GRBs can be identi- 
fied is likely to be a decreasing function of redshift itself. Al- 
though the network of follow-up instruments throughout the 
globe, combined with the aggressive search of high-z GRBs, is 
likely to partly compensate for the difficulties of identifying the 
GRBs at the highest redshifts, such selection effects would still 
have to be carefully folded in to an analysis of the real data (see, 
e.g., Gou et al. 2004 for a recent discussion of the detectability 
of the afterglows of high-redshift GRBs). 

5. CONCLUSIONS 

We find that high-redshift GRB detections are effective at 
constraining the small-scale power spectrum of cosmological 
density fluctuations. Assuming that GRBs trace out the cosmic 
star-formation rate, we generate expected GRB detection rates 
for various cosmological models which include suppression of 
the density fluctuations on small-scales. The effects of such 
suppressions become more notable at high redshifts, where the 
characteristic collapse scales are smaller. Correspondingly, we 
are able to obtain strong constrains from high-redshift GRB 
detections. For example, a Swift detection of a single GRB at 
z ^ 10, or 10 GRBs at z ^ 5, would rule out an exponential sup- 
pression of the power spectrum on scales below R c = 0.09 Mpc 
(exemplified by warm dark matter models with a particle mass 
of m x = 2 keV). Constructing the intrinsic fractional GRB dis- 
tribution from a luminosity-limited sample would provide ad- 
ditional constraints that are independent from the uncertainties 
in the normalization between GRB and star-formation rates, 
and from the uncertainties in the intrinsic luminosity function 
of GRBs. We find that a detection of 1 GRB/yr at z k, 12 also 
provides an upper limit on the running of the spectral index, 
a <; - 0.05, thus potentially placing constraints on the infla- 
tionary potential. 
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and AST-030729 1 (to ZH) and by NASA through grants NAG5- 
26029 and S WIF03 -0000-003 3 (to ZH) and S WIF03 -0020-005 8 
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